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Abstract 

Geodesic acoustic mode (GAM) is analyzed by using modified gyro-kinetic (MGK) equation applicable 
to low-frequency microinstabilities in a rotating axisymmetric plasma. Dispersion relation of GAM in the 
presence of arbitrary Mach number is analytically derived. Toroidal rotation plays the same effects on 
the GAM regardless of the orientation of equilibrium flow. It is shown that the toroidal Mach number M 
increases the GAM frequency and dramatically decreases the Landau damping rate. The valid of classical 
gyro-kinetic (CGK) equation is also examined. For zero electron temperature, CGK is identical with MGK. 
For non-zero electron temperature, CGK gives the same real frequency of GAM as MGK but induces an 
instability with a growth rate proportional to M 3 jq, where q is the safety factor. 
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I. INTRODUCTION 


Geodesic acoustic mode (GAM)[l], seen as the high-frequency branch of zonal flows, is a well 
known phenomenon naturally existing in tokamak plasmas [2]. 


many tokamak devices 


see, for examples, Refs. 


theoretical analyses [6- 13j and numerical simulations 


It was experimentally observed in 
5) and intensively investigated in terms of 


14j-tl7||. GAM is basically a sort of electostatic 


perturbation with toroidally symmetrical and poloidally nearly symmetrical structure. That is, the 
thermal pressure and density perturbations have wave numbers m = 1 and n = 0, while the polodial 
Lagrangian perturbation has a poloidal symmetrical structure with m = 0. The typical frequency 
of GAM in the ideal magnetohydrodynamic (MHD) model reads oj 2 = c 2 (2 + q~ 2 )/R 2 , where q is 
the safety factor and c s = (T(T e + Tj)/mj) 1//2 is the sound speed with the ion mass m;, electron and 
ion temperature T e and Tj, and adiabatic index T. The ideal MHD model yielded T = 5/3 while 
kinetic models predicted uj 2 = (7/4 + r) v\JR 2 QQL where VTi = \/2Ti /mi is the ion thermal 
velocity and r denotes T e /Tj. It was shown that the account of the pressure anisotropy via the 
parallel ion viscosity exactly recovered the adiabatic indexes obtained in kinetic models ( 20 , 21]. 
By using drift kinetic equation, the GAM in the plasma with bi-Maxwellian distribution for ions 
was investigated and the previous kinetic result was recovered when zeroing the anisotropy [22]. 


Recently, Ren studied the GAM in an anisotropic plasma by using the Chew-Goldberger-Low 
(CGL) double entropy equations and found that the kinetic result for zero electron temperature 


was recovered when zeroing the anisotropy 


23]. 


On the other hand, owing to the significant applications to the E x B shear flow control of 
anomalous transport and turbulence, the magnitude, radial profile, and evolution of toroidal flow 
in tokamak plasmas has been an important issue for tokamaks 24^26}]. The equilibrium toroidal 


rotation flow (ETRF) can be on the order of ion thermal velocity 


27H29I] and has been shown to 


El¬ 


be important for GAM and attracted much attention since the seminal monograph by Wang[ 
Wang first investigated the GAM in a toroially rotating tokamak by using the magnetic surface 
averaged poloidal motion equation to eliminate the fast magnetosonic wave. Later, Wahlberg 
found the dispersion relation of GAM by solving the Frieman-Rotenberg eigenvalue equation for 
the Lagrangian perturbation in the presence of ETRF [9] and presented more detailed analysis 
about the GAM and low-frequency MHD modes in a following paper 30j. Lakhin and Ilgisonis used 


the ideal MHD equations in axisymmetric toroidal systems to investigate the GAM by taking into 


account both ETRF and equilibrium poloidal rotation flows (EPRF) 31]. Recently, Ren studied 


the GAM in a toroidally rotating tokamak plasma with an arbitrary (3 by following Wang’s way 
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to deal with the poloidal motion equation [321] and then re-studied the effect of /3 on GAMs in the 


presence of ETRF by considering the harmonics coupling in lieu of magnetic surface averaging 
manipulation J2] . 

The ETRF can increase the frequency of GAM and induce a new low-frequency branch of zonal 
flows 3, |||, and the perturbation is always stable in the MHD model. In view of that the GAM 
has a Landau damping rate 7 oc q 5 e^‘ g2 ^ 4: , it is of great interest to study the collisionless damping 
of GAM in the presence of ETRF, which is the scope of the present work. We also note that 


the GAM driven by the energetic particles can be unstable §[ |l7j, |34j, [35], Considering that the 
plasma rotation flow widely exists in tokamak plasmas, even if no external momentum source is 
injected[36]. In that case, the energetic particles may not start to play a role. Then it is of great 


importance to consider the kinetic effect of ETRF on the GAM. Here with the aid of modified 
gyro-kinetic (MGK) equation { 3 ^, we theoretically analyze the effects of arbitrary toroidal rotation 
on GAM. It is found the Mach number M increases the GAM frequency and remarkably decreases 
the damping rate. The validity of classical GK (CGK) equation is also discussed and shows that 
only for zero r, CGK is identical with MGK. In the case of non-zero r, CGK yields the same 
real frequency of GAM as MGK but introduce an instability with a growth rate proportional to 
r 2 M 3 /q. The rest content is organized as follows. In Section ni the equilibrium distribution 
and MGK equation is presented. The dispersion relation of local GAM in a toroidally symmetric 
tokamak plasma is derived in Section IIIII Section [TV] is devoted to the detailed discussion about 
the GAM frequency and damping rate. The validity of CGK equation is checked in Section [V] and 
finally, the conclusion is performed in Section IVI1 

II. EQUILIBRIUM AND MGK EQUATION 


We consider a large-aspect-ratio tokamak plasma with a toroidally symmetric magnetic field 
B = /(VOVC+VCxVV’, and work in the (r, 0 , () coordinate system, where is the magnetic flux, 
and ( and 6 are the toroidal and poloidal angles, respectively. Circular cross-section is assumed 
in our calculation with R = Rq + r cos 9. The subscript 0 denotes the equilibrium profile and the 
prefix 5 denotes the perturbed one, but the subscript is omitted and the equilibrium magnetic field 
is referred to by B directly. The equilibrium distribution function for ions with an ETRF is given 
by [29 


Fq = n 0 (ip,d)(TTV f. 


2 A-3/2 e -(v-u 0 j 2 /®Ti 


(1) 


3 
















where uq = ujt(iP)R 2 ^C is the toroidal rotational velocity and no is the number density of ions, 


noftM) = iV(^)e M2/(1+r) . 


( 2 ) 


Here, M = ojtR/vti is the Mach number. Noting that T e and T) are both functions of if), M depends 
on the poloidal angle via R. As for electrons, the equilibrium distribution function has a standard 
Maxwellian from, Fq = no(7ru|, e )~ 3 / 2 e~ mev2 A 2Te ). Such a distribution requires r 3> m P Jnm. That 
is, the time scale for electron relaxation can be much shorter than the ion time scale[291]. The two 
distribution functions above are the cornerstones of the analysis. 

The perturbed distribution function is determined by MGK equation, which reads [ 37 , 39. 40] 


Qpi 

5Fj = (8 Fq /dE)qjd(j) + [1 - J o( k rPj)\j^ + Jo^hj, 


( 3 ) 


in which 5hj is governed by 


d_ 

at 


+ (w\\b + u 0 + vd) ■ v 


5hj 


= ~ Qj J o^r(ii + u 0 • v)<^ - 


dE \dt 


rrijui 


-b x V5(j) • VG) 


q ■ p'j 

+J 0 ^rb x V5(j) • [(iO||6 + u 0 ) • Vu 0 + Vu 0 • {w\\b + u 0 )]—- 

Uc 


dE 


( 4 ) 


Only perturbed electrostatic potential is taken into account, which is justified for electrostatic 
GAM in a low-/3 plasma. Here, p = T^rrijw\ is the magnetic moment, E is the energy defined in 
Eq. ([6]) below, k r is the radial wave number of GAM, 5(f) is the perturbed electrostatic potential, 
qj{= ±e) is the charge of species j (j = i,e for ions and electrons, respectively), Jq is the zeroth- 
order Bessel function, pj = wj_/u J c is the Larmor radius, Wc = qjB/nij is the gyro frequency with 
Tfij being the mass, w = v — uq is the particle velocity in the local reference frame moving with the 
velocity uq relative to the lab frame, and vo is the leading order drift velocity. In the lab frame, 
electric field can be expanded as E = E- 1 + Eq + ■ ■ ■, where E—\ = — uq x B and Eq is related to 
d>o. Equilibrium analysis 


29[] shows e<f> 0 = ■ 


Hence in the local reference frame, particles 


only feel the potential <ho- As a result, the drift velocity can be expressed as 


Ed = [(w« + w ±/2)/ UjJ c ]b x VlnB 


b 

+ ~1 x 
LUc 


+ Uq • V-Uq + W\\(b • VU 0 + Uq • V6) 


Meanwhile, the energy E is defined as 


_ t 1 2 1 2 

E = —rrijW - -mjUQ + qj® 0 . 


( 5 ) 


( 6 ) 
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As a result, the equilibrium distribution can be arranged as 


F 3 = N(‘tl>)(‘irv%' i )-V 2 e- E ' T i 


( 7 ) 


III. PERTURBED DISTRIBUTION FUNCTION AND DISPERSION RELATION 


We focus on the ions perturbed distribution function first. Using the properties of Ho, one can 
show that Vuo = ujtR(VRV( — V£Vi?) + i? 2 Vc uq ■ V6 = b ■ Vuoj and 


v D ■'V'lp = + wl/2)V\\ + 


oj^I 


2wj(l + r) 


Vii R 1 


-Vn Rf- 




( 8 ) 


as well as 


(w\\b + u 0 ) ■ Vu 0 + Vn 0 • (w\\b + uq) = 


Iw || 


+ oj t R IVwt- 


( 9 ) 


In view of the two equations above, MGK equation is reduced to 


dgdhi — in l d sin 65h t — i—5hi = —iJ$— — jFqSc/). 

ujI Fiol 


( 10 ) 


t 

Here, = w\\/{qR) is the transit frequency and n l d = k r 5 l b , where 5 l b is the ions orbit width defined 

12 , m2 v 2 


as 


wt. 


+ i} w j_+ 1+ ^‘ +2w\\VTiM ). By assuming 5<f> = the disturbed distribution 


function can be easily solved as 


SFi = --F*[ 1 - J%(k rPi )]6ct> - -J$F* 




ikO 


n,k 


(l - k)6(j) n e 
l — k + uj/ul 


( 11 ) 


Generally for ions, the hnite-Larmor-radius (FLR) effect is taken into account by assuming k rP i ~ 
A < 1 and the finite-orbit-width (FOW) effect is also considered by assuming n d ~ A. While for 
electrons, we have fc r/ o e ~ 0 and n d ~ 0. The electron disturbed distribution function is simplified 
toll 


5F, = 


-F^Wke ike . 


( 12 ) 


k^O 


The dispersion relation of GAM is derived by using quasi-neutrality condition, 5rii = 5n e , 
namely, f d 3 w5Fi = f d?w5F e . Inserting 5Fj into the quasi-neutrality condition will yield the 
dispersion relation. Before doing that, we need to simplify SFj by cutting off the coupling chains. 
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Considering 5<j) n ~ ( k r p l ) n 5(j)Q , we take into account only the coupling between 6<j> o and 5<p±\ by 
neglecting all high-order harmonics. Here, p l = VTi/oj l c is the ion Larrnor radius. Consequently, 
we write 5Fi = 8F® + SF^e^. Eventually, using Eqs. (TTTh and (fT2l) and the quasi-neutrality 
condition, 5cf>k is obtained as [1 + t(1 + CZ(C))\&(t>k = —i\qk r p l TkQ(kM)8^) o for k = ±1. For k = 0, 
we find 


S{M)5<t >0 - g(M)8<k + = 0. 

qk r p l qk r p l 

Finally, we obtain the following dispersion relation 

l g\ M )+g 2 (-M) 

2 r l + r(l + 0Z(C)) ’ 

in which we have denoted 


S(M) = i + | + < 2 + 

2 1 + r 

+ ( 1 + ^ m ' 2 )c 2 + 5 + 

/ i m 2 

G(M) = C + 2M+ U 2 + 2(M + - + — )Z(C). 


m 

c 

M 2 

l + T 


1 + 


M 2 

l+T 


(13) 


(14) 


(15) 

(16) 


Here, ( is defined as quiR/vTi and 2(C) is the plasma dispersion function. The previous result[41] 


is recovered by zeroing M. For convenience of discussion, the Mach number M is assumed to be 
positive since the dispersion relation above is an even function about M. In other words, effects 
of ETRF on the GAM is independent of the fact that the rotational flow is parallel or antiparallel 
to the longitude current as predicted in the MHD model [ 


IV. GAM FREQUENCY AND DAMPING RATE 


Explicit analytical solutions to the dispersion relation (11411 for arbitrary C are difficult to obtain. 
Here we are restricted to the GAM with () 3> 1 to find the asymptotic solution. Since there is 
C 2 = <Z 2 (7/4 + r) in the non-rotating plasma, C ^ 1 requires large safety factor. Hence, analytical 
results below are expected to be accurate enough only when q is high enough. Now we can 
asymptotically expand the plasma dispersion function Z(() = ioyfn rexp (— ( 2 ) — C -1 (l + C~ 2 /2 + 
3C 4 /4 + 15() -6 /8 + •••). Neglecting all terms of order higher than 0(£~ 6 ) leads to the following 
reduced dispersion relation: 


J__0i 

q 2 C 2 


<?0 


c 4 


+ iy/naC 3 e ^ = 0 


(17) 
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in which 


„ 7 r , r2 Af 4 

4 1 + t 

„ 23 9t 3r r , r2 M 4 

g °-¥ + Y + V +5M + 2(TT7) : 


(18) 

(19) 


We now separate ( into two parts, qQx + iQ7d with Qk ^ 7 d- Here, Qk is the normalized 
frequency of GAM and 7 d is the damping rate, f Ik is then determined by 


O 2 - 
il k — 


Gi 


+ 


's+4 


and the imaginary part of (1171) gives the damping rate of GAM in the presence of ETRF as 

V^q 5 ^K 


Id = 




( 20 ) 


( 21 ) 


2^ 2 + 4<5 oA ? 2 

where a = 1 is adopted due to the factor that 7 ^ <C Qk- The two equations above are the major 
results representing the effects of ETRF on the GAM. 

Before further discussion about Eqs. (l20l) and (12TT) . let us pay attention to the case of M = 0. 


The original dispersion relation is reduced to the one in Ref. 


41 


when zeroing M, while the 


simplified dispersion relation (1171) is different with previous results. To illustrate this difference 
clearly, previous results are written here [see, for example, Eq. (7) in Ref. 41 and Eq. (30) in Ref. 
Q as 


c 2 


~o To I T+ r — T7 “^+2r+ — 


c 4 


+ iyfr rC 3 e ^ = 0 . 


( 22 ) 


The coefficients of £ 2 and £ 4 terms are different with their analogs in Eq. (HU). These tiny 
differences are induced by different asymptotic expansion. To obtain a more accurate expression 
of simplified dispersion relation, we multiplied Eq. (fl4l) by 1 + r( 1 + and then kept the 

terms to the leading order. However, on the order of 0(l/q 2 ), Eq. (I17|) yields the GAM frequency 
in non-rotating plasma as 


n 2 a = 


7 


which is identical with the one in Ref. 


+ T 


41 


1 + 


46 + 32t + 8 r 2 
(7 + 4 r) V 


(23) 


Now we take into account the presence of Mach number. Eq. (1201) illustrates the dependence of 
GAM frequency on the safety factor q, temperature ratio r, and the Mach number M. For general 
Mach number M ~ 0(1), by ignoring terms proportional to 1 /q 2 , GAM frequency is reduced to 


o 7 o Af 4 

^ = 4+t + 4M 2 + _-. 


(24) 
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While the MHD model yields (0^ HD ) 2 = T{1+t)+AA1 2 + according to Ref. 9. The coefficients 
of M 2 and M 4 derived in the GK model are identical with their MHD analogs. The only difference 
between the MHD and GK results lies on the adiabatic index T(1 + r) and | + r. Obviously, the 
GAM frequency is increased by the increasing Mach number as shown in the MHD frameworkQ]. 

The toroidal Mach number is shown to decrease the damping rate by increasing the frequency. 
The dependence of the Landau damping rate on the Mach number is plotted in Fig. HI One can see 
that the damping rate is dramatically decreased by the increased M. That is, the Mach number 
tends to increase the GAM frequency and destabilize the GAM by diminishing the damping rate. 
Besides, according to Fig. |TJ the analytical result ()21l) differs from the exact numerical result for 
q = 1 or for q = 2 when M < 0.3. For q > 2 or large M, Eq. (1211) agrees well with the numerical 
result. It is, of course, not surprising since Eq. (1211) is valid for large £, which requires large q or 
large M. 


V. DISPERSION RELATION IN CGK MODEL 


It is argued that CGK model is suitable for describing the toroidally rotating plasmas or not. 

For convenience of discussion, we still focus on ions first. CGK equation gives the ions perturbed 

t2\ dFi 


distribution function as 5F = e(dFQ/dU)5(f> + (1 — Jq)^^ + JoSh with 5h determined by 

(— iu + Vq ■ V)5h = ieJoQ8(f>. 


39|, 


4C] 


(25) 


Here, U = \niiV 2 + e<F is the total energy and Q = uidF^/dU + k x b ■ VFq/^R). The equilibrium 
potential reads <F = d>_i + i n the lab reference frame with 3>_i(^;) = — f UTdip. As a 


rriiOj£,R z 


result, the equilibrium electrostatic field is Eq = cutVt/’ — V • The lowest order velocity is 

Vq = U|| b+ [(vjj + v\/2)/ojl]b x Vln-B + vg. It is easy to obtain that ve ■ W8h = —UTRib ■ V)5h + 
,;i, “t Rs infl x h 

lKr ui(l +T~ 1 ) 0n - 

Poloidal asymmetry of Eq needs to be taken into account due to the definition of U. Noting 
that the spacial dependence of Eq should be calculated at fixed energy U and magnetic moment 
(jl, we rearrange Eq as 


H = N^)e 


u i , 9 , /9 

Ti (nVTi) A 


(26) 


We can find 


**J~2 M^Fj(i + 2l 

Rv I! \ V Ti ^ v Ti 


tM 2 \ 
1 + r ) 


sin@. 


1 + r 


(27) 















M 


FIG. 1: The dependence of Landau damping rate 
on M according to the exact numerical solution to 
the dispersion relation (fT4l) (solid curves) and to the 
analytical expression m, respectively, for r = 1 and 
fixed safety factor, q = 1 (a), q = 2 (b), and q = 3 
(c). 


That is to say, due to the outward shift of ions, the ions distribution function becomes poloidally 
asymmetric under fixed energy U and /r. As a result, we can reexpress Q as Q c — J\f l sin 0 with 
Q c = udF^/dU and 


A fi = 


k r v d Mv Ti 
T i V ll 


K 


(28) 
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The CGK equation (1251) is rewritten as 


de5h — indSmOdh — i — 5h = ieJo—5& — ieJo— smOScj). 

LO t LO t t O t 


(29) 


Here, = k r Vd/cut, Vd = + \v\ — J/(-Ru;*) is the zeroth-order radial drift velocity, 

and oot = (uy /R — cut)/q is the modihed transit frequency. It should be noted that cut = w) and 
Vd = Slcul by recalling uj_ = w± and v» = wn + cutR- After some algebraic manipulation, the 
perturbed distribution function of passing ions is analytically given as 


5Fi 


_JdF& dF* 

~ \dU 

X y1 i n ~ k J n+ i- k {n d )Ji{n d ) 


r) 

l-k 


ikO 


l — k + ou/cut 


+ Ji 


! eJfj 
2icu 


e ie 5(t) - e~ w 5(t> 


-w } 


+2 Y, i"-*-' ( '~ * K " +,' ~ * } J„ + ,- k m n e m 


(30) 


(l ~ k + u/u} t )n d 

For zero M, the total energy U is replaced by the kinetic one E = One then can find 

J\f t = 0 and the perturbed distribution function above is reduced to the previous one in a non¬ 
rotating plasma. As for electrons, there is J\f e = 0. Neglecting the FLR and FOW effects, the 
electron disturbed distribution function is the same as Eq. (El- 
According to Eq. (|26l) . one has 


dFZ_ 
dU " 

M 

Bdfj, 


Ti 
MvTi 
T t v || 


1 - 


Mvn 


U|| 


FL 


(31) 

(32) 


dF e F e 

In view of and after some similar manipulation, we obtain the following dispersion 


relation 


S(M) - ±r 


g(M)V(M) g(-M)V(-M) 


= 0 , 


(33) 


A(M) A(—M) 

in which S{M ) and Q{M) are defined in Eqs. (H5]l and (fl6l) . respectively, and A{M) and V(M) 
are defined as 

( 2 Z(() + M 2 Z(M)\ 


A{M) = 1 + r 1 + 


V(M)=(( + M)(1 + (Z(())- 


C + M 
M 2 




c 


+ 


c 2 Z(() + M 2 Z(M) (1 tM 2 \ 


C(C + M) 


\2 1 + T 




(34) 


(35) 
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Although just like Eq. CHD, this dispersion relation still can reproduce the classical one in the 
non-rotating case, Eqs. (1331) and (1141) are not identical with each other. Only for zero electron 
temperature, CGK is exactly identical with MGK. For simplicity of discussion, let us restricted 
ourselves to the small Mach number case, i.e., the case of M <C 1 since Z(M) is presented in Eq. 
(1551) . Then Z(M) can be asymptotically expanded as iy/i rexp (— M 2 ) — 2M(1 — 2M 2 /3 + •••). 
Thereby, we can suppose the following ordering relation: M ~ e and l/q ~ e. Neglecting all terms 
of order higher than e 4 in Eq. (1331) leads to the following reduced dispersion relation: 

l G'i % 


r 


C 2 


+ iWir 


cr( 3 e 


+ M 3 e ~ M2 T 


(l + 5r- 


2C 2 


c 4 V 2 


(1 + r) 


= 0, 


(36) 


in which 


Q[=‘- + T + ^ + am 2 
4 q 2 


M 4 

+ T77 + 


M 2 /, n2 
— 1 +T 2 
q z 


(37) 


23 r 6 — t 

G ° = ~8 + 4 (1 " 2T) + ~A? T 


+ 


M 2 


(13 - 30r - 15r 2 


(38) 


It should be specifically pointed out that due to coupling effect of 5<f>±\ induced by the nonzero r 
and M, we need multiply the dispersion relation (1331) by £ 2 — M 2 to asymptotically expand the 
equation. As a result, the boxed term in the coefficients above is not accurate, or more precisely, 
can not reduce to the ones in the case of zero t by zeroing t directly. While Eq. (1141) can reproduce 
the dispersion relation with zero t by letting r = 0 directly. 

Albeit there is tiny difference between Eq. (1201) and the real part of Eq. (1361) . the formula above 
yields the real frequency of GAM to the first order as 


O 2 - 

\l k — 


+ r 


1 + 


46 + 32r + 8r 2 


+ AM 2 + 


M 4 

1 + T 


(39) 


(7 + Ar) 2 q 2 

Apparently, this frequency is the same as the MGK result. That is to say, from the aspect of GAM 
frequency, CGK is the same valid and accurate as MGK for small Mach number. When M is on 
the order of unit, the CGK dispersion relation (1361) needs numerical evaluation, which will not be 
presented here. While we pay attention to the damping rate, Eq. (I36|) shows remarkable difference 
with Eq. (1211) . Due to the exponential decay of Landau damping rate, the second term can easily 
exceed the first term in the square brackets of Eq. (1361) . For example, when M > 0.0046 for r = 1 
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and q = 3, the second term becomes greater than the first term. Substituting £ 2 ~ g 2 (7/4 + r), an 
instability comes into being with a growth rate 


7 g 


2-^(1 + 4r)(7 + 4t) 2 


M 3 2 

- T Z 

<1 


(40) 


Due to the fact that 7 ^ in Eq. (12ip and above are both much less than while CGK gives the 
same £Ik as MGK, it can not be analytically determined here which of CGK and MGK is more 
reliable in the case of small Mach number. Experimental observations related to the GAM in the 
presence of torodial rotation may tell the answer. 


VI. CONCLUSION 


Geodesic acoustic mode (GAM) in a toroidally rotating tokamak is investigated by using 
modified gyro-kinetic (MGK) equation. The equilibrium rotational flow is assumed to be along 
the toroidal direction since the zeroth-order poloidal flow is forbidden due to the neoclassical 
constraints[291]. The perturbed distribution function is analytically obtained from the GK equa¬ 


tion with m = n = 0 for the perturbed potential 5(j). By focusing on the passing ions and assuming 
small drift orbit radius, the general dispersion relation of GAM is derived, in which the sonic Mach 
number M represents effects of toroidal rotation and £ is the frequency normalized by VTi/R■ For 
large safety factor, the GAM frequency is f| + r + 4M 2 + j^p to the leading order. It is 
shown that the coefficients of M 2 and M 4 in the kinetic result are the same as the ones in the 
MHD result. Torodial Mach number M increases the GAM frequency as predicted by the MHD 
model[t|, and dramatically decreases the Landau damping rate. The applicability of classical GK 
(CGK) equation for small Mach number, in which Fq = Fq(X, U, /r) and U is defined as -mu 2 + e<f>, 
is also checked. It is found that CGK is identical with MGK for r = 0. In the case of non-zero r, 
only the real frequencies of GAM obtained in both models are the same. The collisionless damping 
term is exceeded and an instability comes into being with a grow rate proportional to r 2 M 3 /g. 
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